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ABSTRACT 

A  fast,  accurate  and  portable  data  reduction  system 
was  developed  for  the  XR-3  Captured  Air  Bubble  testcraft 
being  evaluated  at  the  Naval  Postgraduate  School,  Monterey, 
California.   The  system  consists  of  four  units.   A  magnetic 
tape  recorder  is  used  for  data  play  back.   A  signal  condi- 
tioner unit  with  a  built-in  analog-to-digital  converter  was 
developed  and  is  used  to  filter,  amplify,  sum  and  further 
prepare  the  data  for  transmission  to  either  a  strip  chart 
recorder  or  a  digital  X-Y  plotter  through  a  Monroe  1880 
calculator . 

Preliminary  use  of  curve  fitting  techniques  are  dis- 
cussed; .calculator  programming  and  the  various  problems 
and  solutions  encountered  in  the  development  of  the  system 
are  described. 
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I.    INTRODUCTION 

Associated  with  most  experimental  research  projects 
is  the  requirement  to  obtain,  evaluate,  interpret,  and 
analyze  large  quantities  of  data  in  order  to  draw  meaning- 
ful conclusions.   A  manual  data  reduction  method  is  both 
time  consuming  and  prone  to  error.   On  the  other  hand,  an 
automatic  data  reduction  system  is  more  accurate  and  can 
significantly  shorten  the  man-hours  required  to  reduce 
raw  experimental  data  into  a  form  which  can  be  quickly 
utilized  by  the  project  engineers. 

An  automatic  data  reduction  system  was  needed  for  the 
XR-3  Captured  Air  Bubble  (CAB)  testcraft  being  evaluated 
by  the  Naval  Postgraduate  School  (NPS),  Monterey,  Califor- 
nia.  This  craft  is  being  operated  at  Lake  San  Antonio, 
approximately  100  miles  from  the  school.   Test  data  are 
recorded  on  a  14-track  Pemco,  Model  120B,  magnetic  tape 
recorder.   The  method  used  to  reduce  the  data  has  con- 
sisted of  transporting  the  tape  unit  back  to  NPS,  using 
a  digital  voltmeter  and/or  a  strip  chart  recorder  to  ob- 
tain a  visual  data  presentation,  and  then  making  hand 
performance  plots.   The  ratio  of  data  acquisition  time  to 
reduction  time  was  about  1  to  5.   The  need  for  a  faster 
system  arose  from  the  decision  to  acquire  a  facility  that 
would  permit  the  test  operators  to  remain  overnight  at  the 
test  site.   In  order  to  maximize  the  efficiency  of 


operating  on-site  over  a  several-day  period,  the  results 
of  each  day's  testing  must  be  reduced  and  analyzed  before 
a  schedule  of  tests  can  be  prepared  for  subsequent  test 
runs.   Therefore,  the  time  to  reduce  the  data  must  be 
shortened  significantly. 

The  automatic  data  reduction  system  described  in  this 
paper  fulfills  these  objectives.   It  provides  for  a  better 
utilization  of  operating  time  at  the  test  site  by  provid- 
ing a  quick,  accurate  and  portable  means  of  evaluating  the 
performance  of  the  testcraft. 
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II.   BACKGROUND 

A.   NATURE  OF  THE  PROBLEM 

Both  the  data  acquisition  system  and  the  data  reduc- 
tion methods  for  the  XR-3  have  been  constantly  improved 
since  the  testcraft  arrived  at  the  Postgraduate  School  in 
March,  1970.   Initially,  performance  data  were  limited  to 
thrust,  velocity  and  plenum  cavity  pressure.   These  para- 
meters were  hand  recorded  and  hand  reduced  to  produce 
rough  studies  of  thrust  versus  velocity  and  optimum  plenum 
pressure.   However,  in  March,  1973,  a  new  data  acquisition 
system  was  installed  [Ref.  1]  in  order  to  test  and  evaluate 
the  performance  of  the  XR-3  with  respect  to  various  other 
operating  parameters  such  as  bow  and  stern  seal  pressure, 
turning  speeds,  turning  rates  and  roll  damping.   This  is 
the  present  system  and  the  one  for  which  the  automatic 
data  reduction  system  was  to  be  designed  and  developed.   A 
brief  description  of  the  data  acquisition  system  can  be 
found  in  Appendix  A. 

The  purpose  of  the  new  automatic  data  reduction  system 
was  to  provide  a  quick  and  efficient  means  of  transforming 
raw  tape-recorded  performance  data  from  tests  of  the  XR-3 
into  a  finished  graphical  presentation  capable  of  being 
used  for  evaluating  the  effects  of  varying  operating  para- 
meters.  The  system  was  to  have  the  capability  of  producing 
a  printed  record  of  the  data  converted  into  physical 
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engineering  units,  such  as  knots  for  velocity  vice  the 
millivolts  (raV)  output  from  the  tape  recorder.   The  system 
had  to  be  portable,  in  order  to  be  moved  in  and  out  of  the 
data  reduction  facility,  and  also  small  and  capable  of  be- 
ing operated  in  a  limited  space.   The  equipment  had  to  be 
simple  and  easy  to  use,  requiring  few  inputs  from  the 
operator.   One  of  the  major  considerations  was  that  the 
system  had  to  be  expandable,  which  would  allow  for  the 
implementation  of  future  requirements. 

The  presentation  of  the  results  was  considered  an 
important  factor.   Three-digit  data  accuracy  was  sufficient 
both  for  plotting  and  for  the  printed  output  since  the 
raw  data  are  accurate  only  to  this  order  of  magnitude.   The 
data  would  require  being  scaled  twice,  once  for  the  plot, 
changing  millivolts  to  inches  on  the  graph,  and  once  for 
printing,  if  desired,  changing  millivolts  to  physical  units 
The  system  also  had  to  be  capable  of  performing  operations 
such  as  summing  inputs  from  two  or  more  channels  of  the 
tape  recorder.   One  such  application  is  in  determing  the 
total  drag  on  the  testcraft  at  a  constant  speed.   In  this 
condition  the  drag  equals  the  thrust,  which  in  turn  equals 
the  sum  of  the  thrust  from  both  the  port  and  the  starboard 
engines . 

All  of  the  phenomena  being  tested  occur  at  a  fairly 
low  frequency,  no  greater  than  about  0.5  Hz.   For  this 
reason  the  sampling  rate  need  not  be  much  greater  than 
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one  Hz  according  to  Shannon's  sampling  theorem  (also  known 
as  Kirchhoff's  sampling  theorem).   However,  to  allow  for 
future  expansion  a  faster  sampling  rate  would  be  needed. 

B.   MEETING  THE  OBJECTIVES 

A  Monroe  model  1880  scientific  programmable  printing 
calculator  and  a  Monroe  model  PL4  digital  plotter  were 
chosen  to  be  the  basic  components  of  the  new  data  reduction 
system  (see  Figure  1) .   Because  of  its  small  size  and  rela- 
tively large  instruction  set,  the  calculator  seemed  best 
suited  to  perform  any  manipulations  necessary  to  prepare 
the  data  for  the  plotter  and  also  to  provide  a  paper  tape 
printout  of  the  data  converted  to  the  proper  physical  units 

An  analog-to-digital  (A/D)  converter  was  necessary  to 
transform  the  analog  tape  recorder  signal  to  a  digital  form 
compatible  to  calculator  input.   A  contract  was  let  to 
Santa  Cruz  Engineering  of  Santa  Cruz,  California  in  April, 
1973  to  design  and  build  an  interface  which  would  make  it 
possible  for  the  calculator  to  select  a  tape  recorder 
channel,  then  convert  the  signal  from  that  channel  to  a 
digital  input  to  the  calculator.   Calculator  software  (pro- 
gram) would  convert  this  input  into  the  various  forms 
necessary  to  accomplish  the  remaining  tasks.   A  delivery 
date  of  June  1973  was  set;  however,  an  extension  was  re- 
quired until  late  July  1973. 

As  an  intermediate  method  of  data  reduction  while 
waiting  for  the  interface,  various  methods  of  curve  fitting 
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and  plotting  were  developed  on  both  the  calculator  and  on 
the  IBM  360/67  digital  computer  available  at  NPS.   These 
methods  are  discussed  in  Section  III.   Because  the  calcu- 
lator solutions  were  temporary  in  nature  and  presented 
many  problems,  the  solutions  to  some  of  which  Were  not 
obtained,  only  the  significant  portions  will  be  discussed 
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III.    CURVE  FITTING 

A.   GENERAL  METHODS 

There  are  various  methods  available  to  engineers  for 
fitting  data  to  an  algebraic  relation  that  lends  itself 
to  easy  manipulation  and  plotting  and  that  can  have  a 
physical  significance.   The  methods  explored  for  applica- 
tion to  the  data  reduction  system  were  the  following: 

1.  Piecewise  -  polynomial  interpolation  method 

2.  Difference  method 

3.  Least-squares  approximation  by  orthogonal  poly- 
nomials method 

4.  Least-squares  regression  method 

1 .   Piecewise-Polynomial  Interpolation  Method 

The  cubic  spline  method  of  piecewise-polynomial 
interpolation  was  not  satisfactory  for  use  by  either  the 
calculator  or  the  IBM  360  computer  for  two  reasons. 
First,  the  method  required  knowledge  of  the  slope  of  the 
curve  at  the  initial  and  final  points.   This  information 
is  not  known  and  only  a  gross  estimate  could  be  achieved. 
Secondly,  a  spline  fit  requires  the  curve  to  pass  through 
each  of  the  data  points.   Recording  error  and  signal 
noise  are  not  smoothed  by  this  method  and  as  a  result  an 
inaccurate  model  of  the  physical  event  results.   A  spline 
fit  is  best  suited  to  determine  polynomial  approximations 
for  exact  data  whose  slope  has  discontinuities  [Ref.  2]. 
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2 .   Difference  Method 

The  difference  method  likewise  does  not  lend  it- 
self to  computer  application;  however,  it  is  an  excellent 
tool  for  determining  the  degree  of  fit  required  for  a 
least-squares  type  method.   When  the  difference  method  is 
employed,  the  data  must  be  gathered  such  that  the  depen- 
dent variable  is  obtained  at  equal  increments  of  the  in- 
dependent variable  (X) .   A  difference  table  is  formed 
taking  iterations  until  an  essentially  constant  difference 
is  reached  throughout  a  given  iteration.   The  order  of  the 
polynomial  to  give  the  best  fit  is  the  same  as  the  number 
of  iterations  to  reach  this  constant  difference  [Ref .  3] . 

That  is,  if  the  second  iterations  produce  a  constant  value, 

2 

then  the  model  will  be  of  the  form  y  =  a   +  ax  =  ax 

The  coefficients  of  the  polynomials  are  solved  for  by 
using  the  Gregory-Newton  interpolation  formula: 

r(r-l)(A2y  )    r (r- 1 )  ( r- 2 )  ( A3y  ) 
y  =  yQ  +  r(Ayo)  + 2_  +  _ 


where 


r  = 


x  -  x 


w 


r  (r-1)  ,---,  (r-n)  (A  yQ) 


w  is  the  increment  between  successive  values  of  the 
independent  variable  x,  and 

y  ,  x   are  the  initial  values  of  the  variables 
o    o 
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EXAMPLE  1  DIFFERENCE  METHOD 

As  an  example  of  this  method,  the  data  for  a  third 
order  equation    y  =  5  +  2x    are  presented  in  Table  I. 
Since  the  constant  difference  occurs  in  the  third  itera- 
tion, the  desired  fit  will  be  given  by  a  third  order 
polynomial.   In  addition,  since  the  third  iteration  re- 
sulted in  a  difference  value  with  zero  variance,  the  fit 
will  be  exact. 


TABLE  I 

Successive  Differences  for  Difference  Method  of  Curve 
Fitting . 


Given  Data 
x     y 

Successive  Differences 

1st 
Iteration 
Ay 

2nd 
Iteration 

A2 

A  y 

3rd 
Iteration 

A3y 

4th 
I ter ation 

0        5 
2       21 
4     133 
6     437 
8    1029 
10    2005 

AyQ=   16 

Ayx=  112 

Ay2=  304 

Ay3=  592 

Ay  =  976 
4 

A2y  =   96 

2  ° 

A  y  =  192 

A2y  =  288 

2  2 

Ay  =  384 

A3y  =  96 

3  ° 

A  y  =  96 

A3y2=  96 

aV-  0 
A  y1=  0 

o,  yQ=  5,  w  =  2 


r  = 


x-0 


and 


y 
y 


5  ♦  f(16,+  |(|-l,(|i,+  |(|.1M|.2,(|i,  ♦  0 


=  5  +  2x 
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which  is  the  exact  equation  from  which  the  data  were  de- 
rived.  To  determine  the  order  of  a  model  from  data  which 
is  not  exact,  each  column  of  the  difference  table  must  be 
inspected  for  an  essentially  equal  distribution  of  plus 
and  minus  signs.   The  number  of  iterations  to  that  point 
determines  the  order  of  the  resulting  polynomial  fit. 
The  remaining  procedure  is  as  discussed  in  the  preceding 
example . 

3 .   Least-Squares  Approximation  by  Orthogonal 
Polynomials 

The.  method  of  least-squares  using  orthogonal  poly- 
nomials could  not  be  used  as  a  calculator  solution  to  the 
problem  since  more  storage  is  required  by  the  method  than 
is  available  in  the  calculator.   However,  it  does  provide 
a  very  efficient  and  accurate  solution  when  used  on  the 
IBM  360/67  computer. 

The  method  consists  of  solving  for  successively 
higher  order  orthogonal  polynomials,  such  that  the  sum  of 
which  satisfies  the  requirement  that 


N 


*"*,        n        n 

n=l 


where 


f  (x)  =   exact  solution 
p(x)  =   polynomial  approximation 
N  =   order  of  the  polynomial 
is  as  small  as  possible  [Ref.  2],       A  computer  program  was 
developed  for  use  on  the  IBM  360/67  in  which  the  X  and  y 
data  points,  the  number  of  points  and  the  degree  were 
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entered.   The  output  consisted  of  two  parts:   The  first 
contained: 

1)  The  coefficients  of  the  fitting  polynomial 

2)  An  estimate  of  the  error  of  each  coefficient 

3)  Sum  of  the  square  of  the  deviations 

4)  The  F-ratio 

5)  The  Chi  Square  goodness  of  fit  statistic 

6)  Degrees  of  Freedom  of  the  data,  and 

7)  A  list  of  the  data  points,  the  evaluated  poly- 
nomial at  each  value  of  X,  and  the  differences 
between  the  two  y  values. 

Part  two  is  a  plot  of: 

1)  The  data  points, 

2)  the  fitted  curve,  and 

3)  the  normalized  error. 

The  computer  program  and  sample  output  can  be  found  in 
Appendix  B. 

4 .   Least-Squares  Regression  Method 

The  leas t- squares  regression  method  of  data  fitting 
is  the  most  familiar  and  the  most  widely  used  by  engineers. 
It  is  used  by  hypothesizing  the  existence  of  a  polynomial 
that  will  best  fit  the  data  as  in  the  case  of  orthogonal 
polynomials.   The  degree  of  the  fitting  polynomial  must  be 
determined  by  another  method  or  by  analysis  of  the  results 
in  order  to  determine  the  "best"  fit.   In  the  case  of  an 
Nth  order  fit,  the  fitting  polynomial  has  the  form: 


2       3  N 

y  =  a  +  bx  +  ex   +  dx   +  +  nx 
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Additional  equations  are  needed  to  solve  for  the  coeffi- 
cient a,  b,  c,  d,  ...,  n.  These  are  generated  by  multi- 
plying the  equation  by  x,  x  ,  x  ,  etc.,  then  summing  the 
variable  on  both  sides  to  obtain: 


,y    =  Ma     + 


b  2_j  x    +c  2-»  x    +dLx    + .  .  ,+n  2^j   x 


^xy       =    a  2-»  x    +b^   x         +cL  x         +d^x         +. . ,+n  /]  x 
JLix    y    =    a  2-i  X    +b  2^  x  +c  2-r  x  +d     Y^x  +  .  .  .  +n   /^  x 


n  +  2 


(1) 


En  V*      n    ,    V*      n+1         V1     n+2     ,    V*     n+3  V*     2n 

x    y    =    a^   x    +b/j  x  +c    Z-»  x  +d    z_j  x  +  .  .  . +n    Z^x 


where 


M  is  the  number  of  points. 

In  matrix  notation  Equation  (1)  can  be  written 
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or  as  : 
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EXAMPLE  2   Least-squares  regression  method 

As  an  example,  the  data  of  example  one  will  be  fitted 
to  a  third  order  polynomial.   The  data  in  Table  I  can  be 
used  to  generate  the  data  in  Table  II. 

TABLE  II 
Least-Squares  Data 


X 

y 

2 

X 

xy 

3 

X 

2 
x  y 

4 

X 

3 
X  y 

5 

X 

6 

X 

0 

5 

0 

0 

0 

0 

0 

0 

0 

0 

2 

21 

4 

42 

8 

84 

16 

168 

32 

64' 

4 

133 

16 

532 

64 

2128 

256 

8512 

1024 

4096 

6 

437 

36 

2622 

216 

15732 

1296 

94392 

7776 

46656 

8 

1029 

64 

8232 

512 

65856 

4096 

526848 

32768 

262144 

10 

2005 

100 

20050 

1000 

200500 

10000 

2005000 

10000 

100000 

Z  |  30 

3630 

220 

31478 

1800 

284300 

15664 

2634920 

51600 

412960 

Substituting  the  data  from  Table  II  into  the  Equation  (1) 


3630  =  6a  +  30b  +  220c  +  1800d 

31478  =  30a  +  220b  +  1800c  +  15664d 

284300  =  220a  +  1800b  +  15664c  +  51600d 

2634920  =  1800a  +  15664b  +  51600c  +  412960d 


(4) 


Solving  (4)  for  the  coefficients  gives 

a=5,  b  =  0 ,  c  =  0 ,  d  =  2 ; 
therefore  the  resulting  equation  is 

y  =  5  +  2x 

which  is  the  exact  solution  to  the  data  being  fitted 


22 


B.   CALCULATOR  APPLICATIONS 

Experimentation  with  representative  data  on  the  IBM 
360/67  computer  showed  that  a  polynomial  of  seventh  order 
was  required  to  best-fit  a  curve  smoothly  through  the 
data  points.   This  was  also  verified  using  the  difference 
method.   The  least-squares  regression  method  of  curve  fit- 
ting was  implemented  on  the  calculator  solving  for  a 
seventh  order  polynomial.   Various  problems  were  en- 
countered while  trying  to  solve  Equation  (1)  for  the  un- 
known coefficients  (a,b,c,  ...,  n)  because  the  A  matrix 
was  highly  ill-conditioned.   Two  methods  were  used  to  try 
to  solve  the  equation.   The  first  was  a  direct  approach 
using  triangular  factorization  with  forward  and  back  sub- 
stitution.  The  second  was  an  iterative  method  using  a 
Gauss-Seidel  iterative  scheme  with  and  without  successive 
over-relaxation  (SOR).   Examples  of  these  methods  can  be 
found  in  reference  4. 

The  triangular  factorization  did  not  produce  accept- 
able results  because  of  the  loss  of  significant  digits. 
The  calculator  solution  differed  so  extremely  from  the 
computer  solution  and  the  plot  of  the  resulting  polynomial 
equation  using  the  calculator  solution  bore  so  little 
resemblance  to  actual  performance  plot  that  other  methods 
had  to  be  tried.   Scaling  of  the  input  data  had  no  correct- 
ing effect  on  the  results  using  this  procedure. 

The  Gauss-Seidel  iterative  method  occasionally  would 
provide  the  proper  coefficients  for  the  fitting  polynomial; 


23 


however,  this  was  also  an  unsatisfactory  method  as  a 
general  solution  to  the  problem.   See  Figure  2  for  the 
flow  chart  of  the  calculator  program  implementing  the 
Gauss-Seidel  iterative  method.   The  initial  guess  that 
was  used  to  arrive  at  the  proper  solution  was  extremely 
critical.   An  improper  guess  would  lead  to  instability 
and  the  coefficients  would  diverge  from  the  proper  solu- 
tion.  Even  when  stable  and  converging,  the  time  required 
for  the  calculator   to  reach  a  solution  was  excessive. 
An  iteration  required  about  20  seconds  and,  depending  on 
the  initial  guess,  a  problem  could  require  over  200  itera- 
tions to  arrive  at  a  solution  which  would  closely  approxi- 
mate the  physical  phenomenon.   The  introduction  of  a 
relaxation  factor  helped;  however,  like  the  initial  guess, 
the  choice  was  critical  and  no  way  was  developed  to  aid 
in  the  choice  of  either  an  initial  guess  or  a  relaxation 
factor . 

The  time  problem  and  the  initial  guess  and  relaxation 
factor  problems  caused  this  method  to  be  unacceptable  as 
a  general  solution  to  the  problem  of  curve  fitting  on  the 
Monroe  1880  calculator. 
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SYM  2 


SYM  3 


SYM  6 


SYM  4 


SYM 

CHANGE 
SIGN 


depress  JUMP,  0,  0 
load  cards  1-1  and  1-2 
place  SENSE  SWITCH  up 
depress  JUMP,  SYM,  1,  RESUME 


SYM  1 


enter  X,  depress  RESUME 
enter  Y,  depress  RESUME 


YES 


NO 


enter   initial  guess   for  Gauss   Seidel   iteration 
depress   RESUME 


depress   JUMP,    0,    0 

load  cards   2-1   and   2-2 

place   SENSE   SWITCH  up    (to   iterate) 

depress   JUMP,    SYM,    6,    RESUME 


enter   relaxation   factor 
depress   RESUME 


Figure    2.        Flow    Chart    of    Calculator    Curve    Fitting    Routine 
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© 


NO 


depress  JUMP,  SYM,  9,  RESUME 


enter  X,  depress  RESUME 
calculated  y  is  printed 


YES 


-^ 


depress   JUMP 
SYM,    4,    RESUME 


YES 


depress   JUMP, 
SYM,    6    RESUME 
enter   new 
guess 
depress   RESUME 


Figure    2.        (continued) 
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IV.   DATA  REDUCTION  SYSTEM 

A.    SYSTEM  COMPONENTS 

The  data  reduction  system  is  composed  of  the  following 
five  major  components: 

1.  Signal  selector  and  conditioner  unit 

2.  A/D  converter  and  calculator  interface  module 

3.  Monroe  1880  calculator 

4.  Monroe  PL4  digital  plotter 

5.  Hewlett  Packard  model  7100B  strip  chart  recorder 
The  data  in  the  form  of  an  analog  signal  are  played 

back  from  all  14  channels  of  the  Pemco  magnetic  tape  re- 
corder.  These  signals  go  into  the  signal  selector  and 
conditioner  unit  from  which  the  operator  chooses  the  path 
that  he  wants  each  of  the  signals  to  follow.   He  can  send 
up  to  eight  signals  to  the  A/D  converter  module  also 
located  in  the  signal  selector  and  conditioner  unit,  which 
in  turn  feeds  digital  raw  data  to  the  Monroe  1880  calcula- 
tor for  manipulation  and  plotting,  or  he  can  send  the  raw 
analog  data  to  the  multi-channel  strip  chart  recorder. 
This  flexibility  is  essential  to  the  operation  of  an 
efficient  data  reduction  system. 

1 .   Signal  Selector  and  Conditioner  Unit 

The  signal  selector  and  conditioner  module  was  de- 
signed to  allow  for  operator  ease  in  selecting,  monitoring, 
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and  routing  signals  from  the  tape  recorder  to  the  other 
components  of  the  data  reduction  system.   All  channels  of 
the  recorder  feed  into  the  module  through  the  first  14 
terminals  of  a  CO-ORD  Switch  10  x  15  pin,  2-deck  Matrix 
Program  Board.   By  selecting  the  desired  output  channel, 
the  operator  can  filter,  amplify  and/or  sum  signals  using 
any  of  nine  filter/amplifier  circuits  and  a  summing  cir- 
cuit on  the  10  outputs  of  the  matrix  (see  Figure  3) .   The 
signal  voltage  can  be  monitored  at  any  of  24  points  using 
a  built-in  Datel  model  DM  2000AR  Digital  Panel  Meter  and 
a  24  position  rotary  selector  switch.   The  monitoring 
points  include  each  of  the  inputs  from  the  14  tape  re- 
corder channels,  each  of  the  outputs  from  the  nine  filter/ 
amplifier  circuits  and  the  output  from  the  summing  circuit 

Signal  conditioning  was  needed  to  filter  out  the 
extraneous  high  frequency  noise  from  the  tape  recorder 
signals  which  originated  from  the  XR-3  onboard  sensors. 
The  noise  results  principally  from  vibrations  of  the  test- 
craft  during  normal  operations;  consequently  each  of  the 
signals  has  about  the  same  characteristic  frequency  noise 
as  shown  in  Figure  4.   The  filters  used  are  a  Butterworth 
low  pass  type  [Ref.  5]  with  a  cutoff  frequency  of  0.5  Hz. 
The  electronic  circuit  diagram  and  the  signal  to  frequency 
response  curve  are  shown  in  Figure  5.   The  circuit  was 
constructed  using  a  Philbrick  Nexus  SQ-lOa  amplifier. 

The  output  from  the  filter  is  then  amplified  using 
a  Burr  Brown  3440J  amplifier  circuit.   The  one-volt 
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Figure  3.   Signal  Selector  and  Conditioner  Unit. 
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Figure  4.   Sensor  Output  Signals. 
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Figure  5.   Filter  Frequency  Response  Diagram 
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maximum  signal  from  the  tape  recorder  must  be  amplified 
to  a  five-volt  maximum  signal  for  use  with  the  A/D  con- 
verter.  This  is  to  ensure  the  best  possible  data 
accuracy  when  using  the  A/D  converter  and  calculator  in- 
terface module.   Since  each  of  the  sensor/amplifier/ 
signal  conditioner  circuits  on  the  testcraft  does  not 
maintain  a  null  output  for  a  zero  signal  input,  a  zero 
offset  voltage  circuit  was  included.   Inasmuch  as  a  full 
scale  input  does  not  always  have  an  exactly  1.0  volt  out- 
put, a  variable  gain  resistor  was  placed  in  the  feedback 
circuit  of  the  amplifier  to  allow  the  operator  to  adjust 
the  gain  to  have  a  5.0  volt  output  for  a  full  scale  input. 
See  Figure  6  for  the  amplifier  circuit. 

Two  of  the  filter/amplifier  circuits  are  added 
then  halved  through  a  summing  circuit  using  another  Burr 
Brown  3440J  amplifier  (see  Figure  7).   This  circuit  is 
used  in  obtaining  total  thrust  from  the  individual  port 
and  starboard  thrust  components. 

Electrical  power  is  supplied  to  the  Burr  Brown  and 
Philbrick  amplifiers  from  a  Philbr ick/Nexus  NPS-300  +15v 
-15v  power  supply.   Zero  offset  voltage  is  supplied  by  a 
SRC  MOD  3564  power  supply  and  the  Datel  Digital  Panel 
Meter  is  powered  by  a  Acopian  MOD  10J75  lOv  dc  power  supply 
reduced  to  5v  dc  by  a  National  Semiconductor  LM340K-7805 
5v  dc  regulator. 
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Figure  6.   Amplifier  Circuit. 
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2.   Analog-to  Digital  Converter  and  Calculator 
Interface  Module 

The  analog-to-digital  converter  and  calculator  in- 
terface module,  also  called  the  A/D  interface,  was  pur- 
chased from  Santa  Cruz  Engineering  in  late  July  19730  Many 
modifications  were  required  to  permit  it  to  be  used  with 
the  Monroe  calculator  and  the  other  components  of  the  data 
reduction  system.   These  modifications  are  discussed  in 
Section  V.  , 

The  A/D  interface  is  selected  through  connections 
on  the  front  of  the  signal  selector  unit.   As  many  as  eight 
signals  can  be  connected  and  generally  the  outputs  from  the 
filter/amplifiers  are  used  as  the  inputs.   The  calculator 
program  selects  which  input  signal  is  to  be  converted  by 
sending  a  three-digit  octal  code  from  the  calculator  IX 
register  to  a  TTL/MST  74164,  S-bit,  serial  to  parallel  con- 
verter in  the  interface.   The  first  three  bits  are  then 
sent  to  a  RCA  COS/MOS  analog  multiplexer  CD4051A  for  the 
channel  selection.   The  multiplexer  receives  the  input 
signals  from  all  eight  input  channels.   It  acts  as  a  binary 
controlled  analog  switch  by  selecting  which  channel  is  to 
be  outputted  to  the  A/D  converter.   The  A/D  converter  is 
then  triggered  by  a  TTL/MONOSTABLE  9602  multivibrator. 

The  A/D  converter  is  a  high  speed,  highly  accurate 
and  stable  Analogic  MP  2410  unit  housed  in  a  small  4  inch 
x  2  inch  x  0.4  inch  package.   It  is  capable  of  ten-bit 
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resolution;  however,  only  eight  bits  are  used  for  com- 
patibility with  the  Monroe  1880.   It  is  capable  of  ranges 
of  -10V  to  +  10V,  -5V  to  +5V  and  OV  to  +5V  with  output  in 
one's  complement,  two's  complement,  offset  binary  or 
unipolar  binary.   The  OV  to  +5V  range  using  unipolar 
binary  was  selected  for  the  data  reduction  system  appli- 
cation.  The  output  from  the  A/D  converter  then  goes  into 
a  TTL/MSI  7491,  eight-bit,  serial-in  serial-out,  shift 
register.   Then  upon  another  command  from  the  calculator, 
the  binary  signal  is  transferred  from  the  shift  register 
in  the  interface  to  the  IX  register  of  the  calculator  for 
further  processing. 

Timing  for  the  interface  is  provided  by  the  A/D 
converter  for  the  internal  function  and  by  the  calculator 
for  data  transfer.   Power  is  provided  by  an  Analogic 
AN3001  modular  power  supply  outputting  +15V,  -15V  and 
+5V  dc.   Timing  and  transfer  logic  is  provided  by  FAIR- 
CHILD  TTL/LSI  7400  Quad  two-input  NAND  GATE  and  7410 
triple  three-input  NAND  GATE  integrated  circuits.   The 
logic  diagram  for  the  A/D  interface  can  be  found  in 
Figure  8. 

3 .   Monroe  1880  Calculator 

The  Monroe  Model  1880  Scientific  Programmable 
Printing  Calculator  was  chosen  as  the  heart  of  the  data 
reduction  system  because  it  combines  the  features  of  a 
simple-to-operate  but  powerful  scientific  calculator  with 
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an  extensive  problem-solving  capability.   The  model  being 
used  provides  1024  memory  locations  for  program  steps  (in- 
structions), 128  main  data  registers  and  10  temporary 
(scratch  pad)  registers  (both  for  data  storage).   The  main 
data  registers  can  be  programmed  to  hold  an  additional 
1024  instruction  steps  instead  of  data  if  a  larger  program 
is  necessary.   Data  is  stored  to  an  accuracy  of  13  signi- 
ficant digits . 

The  operations  of  the  calculator  fall  under  six 
general  categories:   input,  output,  storage,  control, 
execution  and  arithmetic  computation.   All  of  these  opera- 
tions can  be  programmed  from  the  keyboard  either  through 
individual  function  keys  or  through  an  ENTER  CODE  key 
which  allows  the  three-digit  octal  code  of  the  instruction 
to  be  entered  into  memory  (see  Figure  9) . 

A  program  can  also  be  stored  using  a  magnetic  card 
device  which  is  built  into  the  calculator  and  enables  pro- 
grams to  be  both  read  into  the  calculator  and  written  on 
to  small  magnetic  cards  each  capable  of  holding  512  program 
steps.   Data  from  the  main  data  registers  can  also  be 
written  on  the  cards  and  conversely  data  from  the  cards 
can  be  read  into  the  main  data  registers.   Data  may  be  in- 
putted to  and  outputted  from  peripheral  devices  through 
the  IX  register  and  an  external  36-pin  connector  to  the 
peripheral.   Programs  can  be  printed  on  paper  tape  with 
the  listing  containing  the  instruction  number,  the  three- 
digit  octal  code  and  a  symbol  to  identify  the  instruction 
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if  one  is  applicable.   Program  steps  can  also  be  inserted, 
deleted  or  changed  without  having  to  re-enter  the  entire 
code.   This  feature  is  a  tremendous  aid  in  correcting 
errors  in  a  program. 

Control  functions  include  such  items  as  sitting 
flags,  printing  data,  and  setting  the  position  of  the 
decimal  point.   Program  execution  is  extremely  easy  with 
the  Monroe  1880.   Symbolic  program  addressing  makes  it 
possible  to  locate  a  program  anywhere  in  program  memory. 
It  permits  branching  or  jumping  to  any  step  in  a  program 
or  to  the  beginning  of  a  subroutine  without  having  to 
know  the  absolute  storage  location  address.   Branching  re- 
turns to  the  point  of  the  branch  after  the  execution  of  a 
reserve;  a  jump  does  not.   Decision  making  is  made  simple 
with  the  use  of  a  sense  switch,  a  programmable  and  a  key- 
board flag  and  through  the  test  of  a  positive,  negative 
or  zero  value  of  the  contents  of  the  main  entry  register. 
Tests  of  these  devices  are  used  to  control  branches  or 
jumps  to  other  parts  of  the  loaded  program. 

The  arithmetic  functions  performed  by  the  calcula- 
tor include  such  standard  keyboard  items  as  add,  subtract, 
multiply,  divide,  square  root  and  reciprocal.   In  addition, 
other  keyboard  items  such  as  sine/cosine,  arc  sine/arc 
cosine,  logarithm,  statistical  summation  and  factorial  are 
included.   Many  more  functions  and  controls  are  available 
through  enterable  instruction  codes  for  which  keys  are  not 
available.   These  include  most  of  the  codes  for  transferring 
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information  to  peripheral  devices,  program  control,  regis- 
ter manipulation  and  other  arithmetic  functions  such  as 
absolute  value  and  arc  tangent.   A  complete  set  of  instruc- 
tion codes  and  operating  instructions  can  be  found  in 
Refs.  6,  7 ,  and  8. 

4 .   Monroe  PL4  Digital  Plotter 

The  Monroe  PL4  Digital  Plotter  is  a  peripheral  de- 
vice used  in  association  with  the  Monroe  1800  series  calcu- 
lators.  Its  purpose  is  to  display  graphically  data  which 
is  passed  to  it  as  a  result  of  program  execution.   The  graph 
origin  may  be  set  anywhere  on  the  plotting  surface  through 
the  use  of  a  calculator  scaling  routine  which  also  formats 
the  data  for  transfer  to  the  plotter.   The  plot  subroutine 
discussed  in  Section  IVB  and  shown  in  Appendix  D,  determines 
whether  the  plotter  is  "busy"  (actually  plotting)  or  is 
ready  to  receive  data.   If  busy,  the  program  loops  until 
the  plotter  is  ready  at  which  time  it  passes  the  pen  command 
and  the  data  for  plotting. 

Instructions  and  data  are  passed  to  the  plotter 
through  the  calculator  IX  register.   The  four  commands 
which  can  be  issued  by  the  calculator  are  presented  in 
Table  III. 

TABLE  III 
Command  Codes  for  Monroe  PL4  Plotter 

CODE  COMMAND 

160  X  data  to  be  passed 

161  Y  data  to  be  passed 

162  PEN  UP 

163  PEN  DOWN 

40 


The  plotter  responds  to  the  PEN  UP  and  PEN  DOWN  commands 

only  when  the  REMOTE  switch  is  depressed  (see  Figure  10). 

In  the  manual  mode  of  operation,  when  the  REMOTE  switch 

is  not  depressed,  the  pen  responds  to  the  commands  issued 

by  the  operator.   These  commands  are: 

WRITE  (pen  down) 

PEN  UP 

MOMENTARY  (spring  loaded  up,  pen  down  while  depressed) 

and  are  issued  through  switches  on  the  plotter  control 

panel.   After  the  X  data  or  the  Y  data  instruction  has  been 

passed  to  the  plotter,  the  calculator  transfers  the  data 

from  memory  to  the  IX  register  for  transmission  to  the 

plotter.   See  Table  IV  for  further  plotter  characteristics. 

TABLE  IV 

Monroe  PL4  Digital  Plotter  Characteristics 

Plotting  Area  10  inches  x  10  inches 

Ballpoint  Pen  Width          0.009  inch 

Resolution  (X  and  Y  axis)   1000  counts 

Repeatability  0.1  percent  of  full  scale 
Accuracy 

Static  0.2  percent  of  full  scale 

Dynamic  3  inch  retrace  will  overlay 

trace  within  0.035  inches 
between  center  lines 
Dimens  ions 

Width  17  inches 

Length  21  inches 

Height  7.75  inches 

Weight  48  pounds 

5 .   Hewlett  Packard  Model  7100B  Strip  Chart  Recorder 
The  Hewlett  Packard  7100B  Strip  Chart  Recorder 
(Figure  11)  has  two  independent  pen  drives  and  input  modules 
which  can  accommodate  analog  input  up  to  100  V.   The  chart 
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transport  system  has  12  selectable  speeds  and  uses  a 
standard  chart  roll  with  a  10-inch  writing  width.   The 
pens  move  horizontally  across  the  top  of  the  chart  each 
capable  of  a  full  10-inch  travel.   Specifications  are 
listed  in  Table  V. 


TABLE  V 


HP  7100B  Strip  Chart  Recorder  Specifications 


Response  time 
Chart  Speeds 


Input  Range 
Input  Resistance 
Dimensions 

Width 

Height 

Depth 

Weight 


0.5  seconds  full  scale 
1.2  in/h;  0.1,  0.2,  0.5, 

2  in/min 

0.1,0.2,0.5,1,2  in/s 
Variable  from  lmV  to  100  V 
1  megohm 

16-3/4  in 
8-11/16  in 
7-1/4  in 
20  pounds 


B.   CALCULATOR  PROGRAMMING  AND  DATA  TRANSFER 
1 .   Memory 

The  calculator  memory  is  divided  into  four  pages. 
Each  page  is  then  further  divided  into  columns  with  each 
containing  32  registers,  numbered  in  octal  from  00-37. 
The  first  page,  Page  0,  is  user's  programmable  memory  and 
contains  four  columns.   Each  column  can  hold  256  instruc- 
tion codes  giving  a  total  capacity  of  1024  instruction 
storage  locations.   The  second  page  (Page  1)  is  data 
register  memory  also  containing  four  columns  for  a  total 
of  128  data  registers  used  for  main  data  storage.   The 
third  and  fourth  pages  contain  16  columns  each  consisting 
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of  working  storage  and  read-only  memory  (ROM).   The  ROM 
contains  the  permanently  stored  programs  that  perform  the 
keyboard  functions.   The  last  two  columns  of  the  fourth 
page  (Page  3)  containing  the  working  storage,  and  all  of 
the  first  two  pages  are  made  up  of  random  access  memory 
(RAM)  which  can  be  written  into  and  read  from.   Informa- 
tion stored  in  Pages  0,  2  and  3  is  treated  as  a  program 
instruction.   A  single  instruction  (or  program  step) 
occupies  one  byte  (8  bits).   Eight  single  instructions  can 
be  stored  one  item  per  register.   The  format  is  shown  in 
Figure  12.   Each  byte  contains  two  Binary  Coded  Decimal 
(BCD)  digits  except  the  sign  position. 

The  calculator  contains  three  programmable  control 
registers,  the  Entry  (E)  Register,  the  data  Address  Regis- 
ter (DC)  and  the  Index  (IX)  Register.   The  entry  register 
is  the  primary  data  register.   Inputs  from  the  keyboard  are 
put  in  the  E  register.   Output  data  are  formatted  there, 
then  transferred  to  a  print  buffer  and  at  the  end  of  any 
calculation  the  E  register  contains  the  result.   The  entry 
register  is  formatted  the  same  as  a  data  register  with  the 
addition  of  an  overflow  and  an  underflow  digit.   The  address 
(DC)  register  contains  14  bits  and  contains  the  address  of 
a  specific  location  in  memory.   It  is  used  with  various 
macro  instructions  for  defining  the  address  to  be  used  for 
data  or  code  manipulation. 

The  index  (IX)  register  is  a  general  purpose  8  bit 
register.   It  is  divided  into  two  4  bit  parts,  IH  (index 
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high)  and  IL  (index  low).   The  IX  is  used  as  a  buffer 
register  for  the  input  and  output  (I/O)  of  information  to 
peripherals  through  the  I/O  bus.   All  information  trans- 
fer between  peripherals  and  the  calculator  occurs  through 
the  IX  register. 
2  .   I/O  Lines 

The  data  output  line  from  the  calculator  to  the 
I/O  bus  is  BYOT  (byte,  output,  pin  25  on  the  connector). 
Data  from  the  IX  register  is  transmitted  on  this  line  on 
the  falling  edges  of  the  clock  pulses  from  T8-T15  time 
(see  Figure  13)  during  an  OPCB  (output  control  byte,  code 

314)  command  and  during  an  OPDB  (output  data  byte,  code 

315)  command.   The  OPCB  command  is  used  to  transmit  selec- 
tion and  control  data  of  peripheral  devices  while  the 
OPDB  command  will  send  numeric  data  over  the  output  line. 

The  data  input  line  into  the  calculator  from  peri- 
pherals is  BYIN  (byte,  input,  pin  26).   Data  is .transmitted 
from  the  peripheral  when  an  IPDB  (input  data  byte,  code 
305)  command  is  executed. 

The  calculator  transmits  two  signals  to  the  peri- 
pheral device  to  allow  the  device  to  determine  which  type 
of  command,  control  or  data  I/O,  the  calculator  is  execut- 
ing.  The  signal  OTCT  (output  control  time,  pine  23)  re- 


mains at  a  logical  0  from  T8-T15  time  and  DTEN  (data 
transfer  enable,  pin  24)  goes  to  a  logical  1  when  the 
calculator  executes  an  OPCB  (output  control  byte).   If  an 
OPDB  (output  data  byte)  or  an  IPDB  (input  data  byte) 
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command  is  executed,  DTEN  will  remain  at  a  logical  0  from 
T8-T15  time  and  OTCT  will  go  to  logical  1. 
3 .   A/D  Interface  Data  Transfer 

The  operator  determines  which  of  the  input  signals 
to  the  A/D  interface  are  converted  and  sent  to  the  calcu- 
lator by  entering  the  proper  code  (shown  in  Table  VI)  into 
the  computer  program  at  decimal  step  numbers  130  and  135 
(see  Appendix  D).   During  execution,  the  calculator  loads 
the  IX  register  with  the  code  for  the  channel  of  the 
dependent  (y)  variable  (decimal  steps  129  and  130)  then 
outputs  a  control  byte  (OPCB  code  314)  to  set  the  multi- 
plexer to  the  proper  channel  then  convert  the  signal.   The 
next  step  (IPDB  code  305)  transfers  the  converted  signal 
to  the  IX  register.   These  data  are  then  stored  and  the 
code  for  the  channel  of  the  independent  (X)  variable  is 
loaded  into  the  IX  register  and  the  process  is  repeated. 
After  the  two  signals  are  in  the  calculator,  they  are  con- 
verted from  binary  to  BCD,  scaled,  and  printed  if  desired, 
then  passed  to  the  X-Y  plotter. 

TABLE  VI 
A/D  Interface  Channel  Codes 

CHANNEL  CODE 

1  340 

2  140 

3  240 

4  040 

5  300 

6  100 

7  200 

8  000 
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4 .   Digital  Plotter  Data  Transfer 

The  digital  plotter  has  a  slightly  more  compli- 
cated data  transfer  scheme  because  of  its  slower  response 
time.   The  first  42  instructions  (000-04L)  comprise  the 
plotter  subroutine  and  are  used  to  control  data  input  to 
this  device.   The  main  program  branches  into  the  sub- 
routine at  the  decimal  steps  000,  004,  008  and  012.   When 
a  value  is  to  be  passed,  rather  than  a  pen  command,  the 
value  is  stored  in  the  E  register  and  is  register  27  of 
the  last  column.   For  example,  when  transfer  of  Y  data  to 
the  plotter  is  desired,  the  program  branches  to  step  004. 
The  program  in  Appendix  E  can  be  followed. 

The  code  for  Y  data  (161)  is  loaded  into  the  IX 
register  than  the  subroutine  jumps  to  octal  step  16.   The 
contents  of  the  IX  register  is  then  copied  into  another 
register.   The  IX  retains  a  copy,  then  the  calculator  out- 
puts a  control  byte  (OPCB) .   The  program  then  checks  the 
lower  4-bits  of  the  IX  register  (IL)  to  be  not  equal  to  0. 
Since  the  IL  register  contains  a  1,  a  jump  to  octal 
address  25  (decimal  address  21)  is  performed.   A  RINS 
(read  input  status)  is  executed  and  the  plotter  returned 
a  1  to  the  H8      (left  most)  bit  of  IX  if  busy  and  the 
calculator  repeats  until  not  busy  or  the  plotter  returns  a 
0.   The  IX  register  is  then  re-loaded  with  the  initial  Y 
code  (161)  destroyed  by  the  last  step  then  tested  for  a 
pen  command.   This  test  fails  so  the  sign  digit  of  the  E 
register  is  passed  to  IL  then  transferred  to  the  plotter 
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(OPDB) .   The  address  register  (DC)  is  loaded  with  the 
address  of  register  27  of  the  last  column  by  the  instruc- 
tions 273  and  277.  This  register  contains  the  number  in 
the  E  register;  however,  data  transfer  is  easier  from 
this  location.   IX  is  loaded  with  the  contents  of  the 
address  pointed  to  by  DC  which  is  the  two  exponent  digits. 
These  are  transferred  to  the  IX  register  then  DC  is  de- 
cremented by  one  step  (two  digits)  and  the  transfer  is 
repeated.   The  plotter  receives  the  exponent,  the  four 
most  significant  digits  of  the  number  and  the  sign  of  both 
The  plotter  pen  then  moves  to  the  point. 
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V.   A/D  INTERFACE  PROBLEMS  AND  SOLUTIONS 

Initially  the  A/D  interface  had  four  features  which 
inhibited  it  from  being  used  with  the  rest  of  the  data 
reduction  system.   The  two  major  problems  were  (1)  the 
logic  for  data  transfer  and  (2)  the  loss  of  the  most  sig- 
nificant bit  from  the  A/D  converter.   The  two  less  serious 
and  more  easily  corrected  faults  were  (3)  the  voltage  con- 
version range  of  the  A/D  converter  and  (4)  the  conversion 
of  binary  to  decimal  numbers. 

A.   LOGIC  FOR  DATA  TRANSFER 

The  initial  design  of  the  A/D  interface  assumed  both 
DTEN  and  OCTC  to  remain  logical  O's  throughout  the  entire 
normal  time  cycle  instead  of  both  going  high  from  T8-T15 
time.   It  also  assumed  that  DTEN  would  go  high  for  data 
transfer  and  OCTC  would  go  high  for  control  output  instead 
of  missing  the  one  time  pulse  train.   Consequently  the  A/D 
interface  could  not  interpret  the  calculator  requests.   In 
order  to  correct  this  design  flow,  the  circuit  labeled  "A" 
in  Figure  8  was  incorporated.   The  two  signals,  OCTC  and 
DTEN,  are  compared  with  an  exclusive  OR.   If  the  two  are 
different,  the  output  of  the  exclusive  OR  is  high.   The 
original  OCTC  and  DTEN  are  then  "AND"ed  to  determine  if 
either  signal  has  been  missed  indicating  the  command  has 
been  executed.   This  addition  changes  the  OCTC  and  DTEN 
signals  to  that  expected  by  the  rest  of  the  interface. 
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B.  LOSS  OF  MOST  SIGNIFICANT  BIT 

After  correcting  the  first  problem,  analog  signals 
could  be  converted  to  binary  and  transferred  to  the  cal- 
culator.  However,  only  seven  of  the  eight  bits  were 
being  passed  to  the  IX  register.   It  was  determined  that 
the  most  significant  bit  was  not  transferring  into  the 
interface  8-bit  shift  register.   It  was  left  "floating" 
at  the  output  of  the  A/D  converter  then  lost  when  the 
first  transfer  pulse  was  being  issued  by  the  calculator. 
The  circuit  labeled  "B"  in  Figure  8  was  added  to  remove 
the  first  transfer  pulse  during  an  IPDB  and  allows  the 
entire  8  bits  to  be  transferred  into  the  IX  register. 
It  uses  a  TTL/MST  7493  4-Bit  Binary  Counter  and  four  "NOR" 
gates  configured  to  remove  the  T9  time  pulse  from  the 
clock  to  the  8-bit  shift  register  during  the  transfer  of 
data  from  the  A/D  converter. 

C.  VOLTAGE  CONVERSION  RANGE 

The  Analog-to-Digi tal  converter  was  wired  at  delivery 
to  convert  signals  in  a  range  of  -5V  to  +4.96.   Since  all 
the  signals  from  the  XR-3  are  in  a  range  of  OV  to  +1V, 
this  larger  range  reduced  the  significance  which  could  be 
achieved  by  the  converter.   Using  eight  bits  coded  binary, 

the  full  scale  range  of  the  converter  can  be  broken  down 

8 
into  2   or  256  parts.   The  smallest  range  capable  of  the 

A/D  converter  is  OV  to  +4. 98V  so  in  order  to  realize  the 

greatest  possible  accuracy  two  changes  had  to  be  made. 
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First  the  A/D  converter  range  was  changed  to  OV  to  +4. 98V 
and  second  the  signal  condition  and  amplifier  unit  was 
developed  and  built  to  enable  the  tape  recorder  signals 
to  be  amplified  to  a  5V  maximum  signal.   By  making  these 
two  changes  the  full  capability  of  the  A/D  converter  is 
utilized . 

D.   BINARY  TO  DECIMAL  CONVERSION 

The  calculator  data  format  is  binary  coded  decimal 
(BCD)  characters.   Four  bits  are  used  to  represent  a  deci- 
mal digit  and  a  number  is  made  up  of  by  combining  BCD 
characters.   The  A/D  interface  transfers  a  number  to  the 
IX  register  in  binary  representation.   In  order  to  be  used 
by  the  calculator,  this  number  has  to  be  converted  to  BCD 
format.   This  is  done  by  calculator  software  by  a  series 
of  jump  instruction  beginning  at  decimal  step  48  of  the 
data  reduction  program.   IX  register  bit  Ll  (right  most) 
is  the  most  significant  bit.   If  this  bit  is  high,  128  is 
added  to  scratch  pad  register  0.   If  bit  L2  is  on,  64  is 
added,  and  so  on  such  that  if  all  bits  are  on,  scratch  pad 
register  0  totals  255.   After  checking  all  bits  of  the  IX 
register,  the  scratch  pad  register  is  divided  by  256.   The 
result  of  this  is  to  leave  in  the  register  the  decimal 
equivalent  of  the  Ov  to  IV  analog  signal  converted  by  the 
A/C  interface. 
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EXAMPLES  OF  CONVERSION: 

VOLTAGE  IN    A/D  OUTPUT  CALCULATOR  CONVERSION 

0.25V  0100  0000  26/28  =  64/256  =  0.25 

0.50V  1000  0000  27/28  =  128/256  =  0.50 

0.6875V       1011  0000  (  2?  +  25  +  24 ) /28= ( 128  +  3 2  + 16 ) /2 56 

=  0.6875 

1.0V  1111  1111  (27  +  26  +  25+.  .  .  +  2°)/28=  HI-  =.998 
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VI.   CONCLUSIONS 

The  automatic  data  reduction  system  has  been  shown  to 
provide  a  quick  and  accurate  means  of  evaluating  performance 
data  from  the  XR-3  testcraft.   It  has  been  tested  with  both 
simulated  and  actual  tape  recorder  signals  and  proven  to  be 
accurate  to  greater  than  1%  of  full  scale.   Raw  data  can  be 
efficiently  transferred  into  a  finished  graphical  presenta- 
tion and  a  printed  record  of  the  data,  converted  to  physical 
units,  can  be  easily  obtained.   The  system  is  superb  for 
its  simplicity  and  ease  of  operation.   A  job  that  used  to 
require  days  can  now  be  completed  in  a  matter  of  hours. 

The  signal  selector  and  interface  unit  provides  the 
capability  of  zeroing  and  calibrating  the  signals  from  the 
tape  recorder.  It  removes  the  signal  noise  resulting  from 
vibration  and  can  pass  the  data  to  a  strip  chart  recorder, 
or  it  can  internally  digitize  the  information  for  transfer 
to  the  calculator  for  further  processing  and/or  plotting. 

The  provision  for  expansion  has  been  designed  into 
the  system  both  through  the  modification  of  hardware  and 
software.   The  signal  selector  and  conditioner  unit  has 
room  for  internal  modification  and  additions  to  satisfy 
requirements  which  may  arise  in  the  future,  and  the  calcu- 
lator has  over  six  columns  of  unused  memory  which  may  be 
used  to  increase  the  data  reduction  program  to  fulfill 
future  needs  in  this  area. 
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APPENDIX  A 
XR-3  DATA  ACQUISITION  SYSTEM 

The  data  acquisition  system  presently  installed  in  the 
XR-3  testcraft  is  comprised  of  various  sensors  feeding 
through  amplifiers  and/or  signal  conditioners  in  order  to 
prepare  the  signal  for  recording  on  a  Pemco  Model  120-B 
magnetic  tape  recorder  (see  Figure  14). 

A.   SENSOR  SYSTEM 

Sensors  are  presently  installed  to  measure  the  follow- 
ing parameters: 

1.  Port  thrust 

2.  Starboard  thrust 

3.  Forward  seal  pressure 

4.  Aft  seal  pressure 

5.  Plenum  chamber  pressure 

6.  Testcraft  velocity 

7.  Water  emersion  height 

8.  Pitch  angle 

9.  Pitch  rate 

10.  Roll  angle 

11.  Roll  rate 

12.  Yaw  angle 

13 .  Yaw  rate 

14.  Lateral  acceleration 

15.  Longitudinal  acceleration 

Thrust  is  measured  by  a  Revere  USPI-150A  (30)balance 
bridge  transducer.   One  is  installed  on  each  engine  mount 
so  as  to  measure  the  fore  or  aft  thrust  of  each  engine. 


Two  pressure  probes  are  installed  in  the  plenum,  one 
forward  and  one  aft.   Only  one  can  be  used  at  a  time. 
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Figure  14.   Data  Acquisition  System. 
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The  signal  from  the  bridge  is  amplified  by  one  of  ten  Grant 
Model  DCAB-3  amplifiers  to  a  range  between  0.0  and  1.0  volts 
corresponding  to  0.0  to  500.0  pounds  thrust.   Calibration 
signals  of  0.0,  0.5  and  1.0  volts  are  also  available  to 
each  of  the  amplifiers  which  are  found  in  an  amplifier- 
transducer  package. 

Pressures  are  measured  by  pneumatic  probes  connected  to 
Varian  differential  pressure  transducers  accurate  to  +_  1%. 

Pressure  signals  are  amplified  and  calibrated  similar  to 

2 

thrust  so  as  to  have  60  lbf/ft   correspond  to  1.0  volt. 

Velocity  is  measured  by  a  vane  enclosed  in  a  water  probe 
which  transmits  a  frequency  modulated  signal  to  a  Potter 
velocity  conditioner.   The  velocity  conditioner  converts 
the  FM  signal  to  a   dc   signal  from  0.0  to  5.0  volts  corres- 
ponding to  0.0  to  40.0  knots.   This  signal  is  then  reduced 
to  0.0  to  1.0  volt  by  another  amplifier  in  the  amplifier- 
transducer  package.   Calibration  is  performed  by  operation 
over  a  measured  course. 

A  Humphreys  Model  CF18-0101-1  gyro  package  measures 
pitch  angle  and  rate,  roll  angle  and  rate,  yaw  angle  and 
rate,  and  lateral  and  longitudinal  acceleration.   Presently 
all  modes  except  acceleration  are  being  used.   Outputs  are 
between  0.0  and  1.0  volt  and  range  is  +_    20  degrees  in  pitch, 
+  15  degrees  in  roll,  +_    180  degrees  in  yaw,  and  +_  30  degrees 
/sec  in  all  rates. 
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Height  is  measured  by  a  Western  Marine  Electronics 
Model  LM4001  height  sensor  with  the  transducer  located 
forward  of  the  bow.   Range  is  selectable  +_   2    feet  or 
+^  3  feet,  both  with  an  accurate  of  1%.   The  +_  2    feet 
range  is  presently  in  use.   The  signal  from  the  sensor 
is  converted  to  0.0  to  1.0  volt  dc  by  a  height  sensor 
conditioner. 

B.   DATA  RECORDER 

The  conditioned  signals  from  the  various  on  board 
sensors  are  sent  to  and  recorded  on  a  Pemco  Model  120-B 
magnetic  reel-to-reel  tape  recorder.   The  recorder  has  14 
channels  for  electronic  data  plus  an  edge  track  for  audio 
recording.   The  unit  can  be  controlled  from  either  the 
pilot's  instrument  panel  or  from  the  unit  itself,  located 
just  aft  of  the  pilot's  cockpit.   The  recorder  is  portable 
weighing  only  100  pounds  and  measuring  9-1/2"  x  17-3/4" 
x  26-1/2".   It  is  easily  installed  and  removed  from  the 
testcraf t . 

The  recorder  can  operate  at  tape  speeds  of  60,  30,  15, 
7-1/2,  3-3/4,  and  1-7/8  inches  per  second  and  uses 
standard  or  precision  NAB  10-1/2  inch  reels  and  1.0  or 
1.5  mil  one-inch  tape. 
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APPENDIX  B 
IBM  360/67  COMPUTER  DATA  FITTING 

A  fortran  program  was  developed  for  use  on  the  IBM 
360/67  computer  to  aid  in  calculator  software  design. 
The  program  consists  of  five  basic  parts  as  follows: 

Part  1:   Initialization  and  control 

Part  2:   Data  point  plotting 

Part  3:   Determination  of  fitting  polynomial  coeffi- 
cients and  error  criteria 

Part  4:   Error  normalization  and  plotting 
Part  5:   Evaluation  and  plotting  of  final  curve 
The  plotting  subroutine,  DRAW,  is  used  to  provide  the 
final  graphical  output.   It  is  available  as  a  callable 
subroutine  and  utilizes  a  Calcomp  Plotter. 

The  subroutine  LSF  provides  the  fitting  polynomial 
coefficients  and  all  the  written  output.   It  is  a  modifi- 
cation of  the  subroutine  LSQF2  and  the  procedures  and 
algorithm  for  the  method  can  be  found  in  Ref.  2. 
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c 

C  PURPOSE 

C  TC    FIT     A    CURVE    TO    A    SET    OF    DATA    PCINTS    USING    THE 

C  LEAST-SCUARES    METHCC    WITH    ORTHOGONAL    POLYNOMIALS 

C 

C  INPUT 

C      1,  N=NUMBER  CF  PCINTS  TO  FIT  (FCRNAT  12) 

C      2.  THE  DATA  UP  TO  60  PCINTS  (FORMAT  2F10.3) 

C  3.  THE  CEGREE  OF  FIT  (1  THRU  20  FORMAT  12) 

C  4.  TITLE  FOR  THE  OUTPUT  AND  PLOT  (FI^ST  48  COLUMNS  CF 

C         2  CARDS;  FIRST  CARD:  TITLE,  SECOND:  NAME) 

C 

C       4  +********  4****=:**  ********  *:$.**>***  ********  **?*****  4*4** 

C 

c 

OIMENSI CN  FPTC2  00) ,XPT(200) ,X1 ( 30 ,2) , Fl ( 30  ,  2  )  ,X2 ( 60  ) 
X  ,CEL(6C) 

REAL*8  TITLE( 10) , LABEL/1         • / ,TI TLE2 ( 12) ,DMAX 

REAL* 8  X(60), F(6G),Y(60),WI(60) ,DELY(60) ,E(21),SB(21) 

DATA    WI/6C*1.0D0/» LABEL/1     E RR« /  ,  LAEE2/ •  •/ 

CATA    TITLE/10*'  •/ tTITLE 2/12*"  •/ 

LAST=0 

K  =  5 

L  =  6 

READ    (K,91)N 
C 

C    TC  READ  IN  THE  DATA  POINTS 
C 

READ(K,93) (X( I) ,F(I) ,I=1,N) 

DC  7  1=1,30 

X1(I,1)  =X(I) 

F1(I,1)  =  F(  I) 

X1(I,2)-X(  1+30) 

F1(I,2)=F(  I  +  3C) 

X2(I)=X(I) 

X2(I+30)=X(I+30) 

7  CONTINUE 

20    READ  (K,91)KM 

IF  (KM.GT..25JG0  TC  30 
KM=-KM 

8  1  =  1 

READ  (K,200)TITLE2 
DC  12  J=l,6 
TITLE( J)=7ITLE2(J) 
12    CONTINUE 

9  NPTS=30 

IF  (N.LT.NPTS)NPTS=N 

IF  ((I. EQ.2) .AND. (N.LT.60) )NPTS=N-30 

CALL  CRAW(NPTS,X1(1,I) ,F1(1,1)  ,1 ,  1 , LABE L  ,  TI TLE2 , 
X.2,  -2, 1,0, 2, 2, 5, 6,0, LA ST) 

IF  (N.LT. 30)1=2 

1=1  +  1 

IF     (I.LT.3JG0    TO    9 
C 
C       TC    PLCT    THE    ERROR 

CALL    LSF( NfKM,XfF,WltYfDELYi B, SB, TITLE) 

DMAX=0. C 

CC  36  1=1, N 

IF  (DABS(CELYU)  )  .GT.OMAXJ  DMAX  =  DAES  (DELY  {  1)  ) 

36  CONTINUE 

DC  37  1=1, N 

DEL  (I) =DELY(Ii/DMAX/10. 

37  CONTINUE 

CALL  DR/>W(N,X2,DEL,2,0,LABEl,TITLE2t 
X.2,.2,1,0,2,2,5,6,C,LAST) 


FIGURE  15c  FORTRAN  CURVE  FITTING  ROUTINE 
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C    TC  PLCT  THE  FUNCTION  CURVE 
C 

XNT  =  (X(M-X(l))/200.0 
XPT(1)=X(1) 

FPT(l)=POLY(KM,XPT,B, 1) 
DC    10    1=2,200 
XPT(I)  =  XPTU-1)+XNT 
FFT(I)  =  PCLY(KMtXPT,6,  I) 
10         CCNTINUE 

CALL    DRAW(200,XPTfFPTt3f0fLABE2 ,TI TLE2 ,EXSC  , 
X. 2,. 2, 1,0, 2, 2, 5, 6,0, LAST) 
GC    TO    20 

90  FORMAT  (LF10.4) 

91  FCRMAT  (12) 

92  FORMAT  (13) 

93  FORMAT  (  2 FJ. 0.3) 

98    FORMAT  (•  CMAX=« ,F10.4) 
200   FORMAT  (6A8) 
30    STOP 
ENC 


C     

c 

C  FUNCTION  POLY 

C  PURPOSE 

C  TC  EVALUATE  THE  POLYNOMIAL  USING  NESTING 

C  INPUTS 

C  1.    N=THE    DEGREE    OF    THE    PCLYNCMIAL    FIT 

C  2.    X=    THE    ARRAY    OF    POINTS    AT    wHICF    TO    EVAL    THE 

C  POLYNOMIAL 

C  3-  B=  THE  ARRAY  OF  COEFFICIENTS 

C  4..     I=THE    SUBSCRIPT    OF    THE    X 


C 


FUNCTION  POLY(N,X,B,I) 

REAL=?8  6(21),  DX,  C 

DIMENSION  X(200) 

CX=X( I) 

C=6(N+1) 

DC  5  K=1,N 

J=N-K+1 

C=6(J)+C*DX 

CONTINUE 

PCLY=C 

RETURN 

END 


FIGURE  15o(C0NTINUEC) 
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c 

c 

C  SUeROUTINE  LSF 

C  PUPPCSE 

C      70  SOLVE  FOR  THE  LEAST-SQUARES  POLYNOMIAL 

C  INPUTS 

C      1.  ^=NUPBER  OF  OATA  PCINTS(LESS  TFEN  100) 

C      2.  KM=TFE  DEGREE  OF  FIT 

C  3.  X=ARRAY  OF  X  VALVES 

C  4.    F2=ARRAY    OF    Y    VALUES 

C  5.    WI=ARPAY    OF    WEIGHS    FOR    THE    DATA    POINTS 

C  6.    Y=GUTPUT    ARRAY    CF    The    FIT    DATA 

C  7.    DELY=CUTPUT    ARRAY    OF    ERRCR     (F2-Y) 

C  8.    B=OUTPUT    ARRAR    OF    THE    COEFFICIENTS    OF     THE 

C  POLYNOMIAL    FIT 

C  9.    SB=OUTPUT    ARRAY    OF     THE    ESTIMATED    ERROR    IN    B 

C  10.    TITLE=FEACING    FOR    THE    OUTPUT 

C  NOTE:     XtF2fWIfYiDELYfBtSBfAND    TITLE    ARE    REAL*8 

C 

c 

SLBROUTINE    LSF( M, KM t X,F2, WI ,Y,DELY,Bt SB, TITLE) 

REALMS    X(l),F2(l)tWI(l)tY(l)tDELY(l)tB( 1AV SBll) t 
X         FMfFNR*PXF  ,PXP.XPXPtXPXPM,PPXPPfXP,ALPFAfBtTAfPPXF, 
X    F(100)  ,P(  LOO)tP.MI  ICO)  ,S(  LOO)  ,ST(  ICO)  ,A 12  1,21)  ,T(10C)  , 
X    SIC-2,SIG3,SUM£V2,  FLEV  ,  F2  3AR,  SQ  ,  FMF  ,  FMKF  ,  £K  rCHl ,  *  (  ICC) 
X     .TITLEUI  tF3AR,XB-AR 

CC    1    1=1,21 

DC    1    J=l,21 
1    A(ItJ)=C.O 

A(l,l)=1.0 

A(2,2)=1.0 

PBAR=O.C 

SUMEV2=C. DO 

FN=0.0 

DC    5    I  =  1,M 

IF(WI(I  ).EC.0.0)    WI(I)  =  1.0DC 
5    FP*FM*WIU) 

FMR=1.0D0/FM 

FE£R=0.0 

XEAR=0.0 

DC  10  I  =  1,M 

W(I)=WI  (I)^FMR 

PMI)=DSGRT(,M(I) ) 

FU)=F2U  )>PM(I) 

SUMEV2=SU-1EV2+F2(  I)+WHI) 

FBAP=FBAR+F(I  )-*PM{  I) 
1C    XBAR=XSAR+X<I)*WUJ 

T  (1)  =  FBAK 

A(2,1)=-XBAR 

PXF=0.0 

PXP=0.0 

F2BAR=SUMEV2/FM 

SUMEV2=0.D0 

DC  20  1=1 , N 

SLMEV2=SUMhV2+  W(  I  )+ ( F2 ( I ) -F2BAP )**2 

P(I)=(X(I)-XBAR)*PM(I) 

P>F=PXF+P(  I)^F(I) 
2C  P>P=PXP+P( I)*P< I ) 

fFP    =    F*    +    i.C-5 

T(2)=PXF/PXP 

Pr<XPM=l  .  GOC 

s  d)  =  py,xPM 

KNP  =  IASS(K//)  +  l 

B(1)  =  T(  1)-A(1,1)+T(2)*A(2, 1) 

B(2)=T(2)*A(2  ,2) 

DO    1000    K=2,KMP 

CH=0. 

KM=K-1 

Kf2=K-2 

FIGURE  15o(C0NTINUED) 
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If  (KM2.  EQ„0)    GO    TC    200 
XPXP=O.C 
XPXPM=0.0 
B(K)=O.C 
CO    50    J  =  1,M 
XP=X( J)*P< J) 
XPXP=XPXP+XP*P( J) 
50    XFXPM  =  XPX.PM+XP*PM(J) 
ALPHA=XPXP/PXP 
BETA=XPXPM/PMXPM 
PFXF=0.0 
FPXPP=0.0 
DO    100    1=1, M 

pt=pu) 

PU)  =  (X(I )-ALPHA)*P(I)-BETA*PM( I) 

PPXF=PPXF+P(I)*F(I) 

PPXPP=PPXPF+P(1 )*P(I) 
100  PMI)=PT 

T(K)=PPXF/FPXPP 

PMXPM=PXP 

FXP=FPXPP 

A(K,1 )=-ALPHA*A(KMl ,1 } -BETA* A ( KV2 , 1 ) 

A(KfKMl)=A(KMltKM2>-A(KMl,KMl)*ALPHA 

A(KfK)  =  l.jQ 

IF(K.LE.3)  GO  TO  150 

CC  12J  1=2, KM2 
12 C  A(K,I)=A(KM1  f  1-1  J-AlPHA*A(KM1,  I  )-EETA=r-A  (KM 2  ,  I ) 
15C  DO  160  1=1, K 
160  E(I)=B(  XJ.+T(K)*A(Kf  I) 
2CC  SIG3=0.0 

DC    220    1=1,  M 

SC=E(K) 

KKC=K-1 

DO    230    IQ=1,KKQ 

KMC=K-IQ 
23C    SC  =  XU)*SQ+B(KMIC) 

Y(I)=SQ 

C£LY(I)=Y(I)-F2(  I  ) 
22C    SIG3=SIG3-*:    W< I )*OELY( I )**2 

FM<F=MFM-K 

SIG2=SIG3*FM/FMKF 

FLEV=( SCMEV2-SIG3)/SIG2 

SLMEV2=SIG3 

S(K)=PXP 

CC    240     1=1, K 
24C    ST(I)=SIG2/S(I) 

DO    300    1  =  1, K 

SE(I)=G.C 

OC    2p0    J=1,K 
250    SE(  I)  =  5B(I )+ST(J)*A(J,I)**2 
300    Sem=DSQRT(SB(in 

IF(KM.GT.G)GO    TO    301 

CONT  INUE 

IF(K.LT.KMP)GO    TO     1000 
3C1    CONTINUE 

WRITE(o,9510J       (TITLE(  I), 1  =  1, 10) 
9510    FCRMAT(1H1,////,30X,10A8,// ,25X, 'COEFFICIENTS    OF     THE', 
X'    POWER    SERIES    E XPANSICN1 , //,26X, • YCXI =B{ 1)*B (2 )*X« , 
X,+B(3)*X*'*2+B{4)*X**3+...1//) 

WRITE(6,9520)     (  < I , 8 ( I ) ) ,  I  =  1  ,K ) 
9520    FCRMAT((15X,3(«       B I  '  , 12 , • ) = • , 1PD12 . 5) ) ) 

WFITE(6,9530) 
9530    FCRMAK //,25X,  'ESTIMATES    CF    ERRCR    FOR    THE', 
X'     COEFFICIENTS',//) 

WRITE (6  ,9540) (( I,SE(I) )  ,i  =  l,K) 
954C    FORMAT!  (15X,3( '        ERRB ( '  ,  I  2  ,  ' ) ='  ,1 PC9  .3 )  )  ) 

WFITE(6,9570)    SIG3,FLEV 
9570    FCRMAT( //,21X,' SUM    SQ    OEV    = • , 1PD1C. 3, 11X , ■ F-RATI 0    =' 


XtlPD10.3) 


FIGURE    15.  (CONTINUEDJ 
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00  25  IRI=1,M 
25  CFI=CHI*0ELYURI)*DELY(IRI)*WI(  IR  I  l/DABS  (Y(  IRI  i  ) 

WRITE (6,1111)CHI tFMKF 
1111  FORMAT!  «C« ,20X, •CHISU=' tlPD11.3  ,15X,*DEG  CF  FREEDC^  =■ 
>,CPF4.0) 

WRITt(6,955G) 

9550    FCRMAT(//,34X,' DATA FIT  DIFFERENCE1,/ 

2  30X,«I  X(I)  Yd)  YY(I)  YYIIJ-Yd)',//) 

WRITE (6,9  560) (( I ,X( I J ,F2(I)  ,Y(I)  ,CELY(I) )  ,I=1,M) 
9560    FCRMAT(15Xf I16f3F8.3f F11.3) 
10CC    CCNTINUE 

WRITE(6,9580) 
9580    FCRMATC'L1 ) 

RETURN 

ENC 


FIGURE    15. (CONTINUED) 
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VELOCITY    VS    THRUST    5TH    ORDER    FIT 

COEFFICIENTS    OF    THE    POWER    SERIES    EXPANSION 

Y(X)=B(l)*e(2)*X+B(3)*X**2+B(4)*X**3+... 


e(  1J=  2.57974D-03   B( 
B(  4)=  3.2S749D  02   B( 


2)=  1.101340  CI 
5)=-5.34132D  02 


B(  3)=-9.198S6D  01 
B(  6J=  3.26C51C  02 


ESTIMATES  OF  ERRCR  FCR  THE  COEFFICIENTS 


ERRB(  1)=8. 7260-02 
ERRB(  4)=2.816D  C2 


ERRB(  2)=4.646D  00 
ERRB(  5)=5.556D  02 


ERRB(  3)=6.CC4D  01 
ERRB(  6)=3.913C  02 


SUM  SQ  DEV  =  3.0510-04 
CHISQ=   1.271D-C2 


F-RATIG    =    6 .944C-01 
DEG    OF    FREEDOM    =      9. 


. 

— i — DATA 

FIT 

CIFFERENCE 

I 

X(I) 

Y(I) 

YY(I) 

YY(I)-YU) 

1 

0.0 

O.C 

0.003 

0.003 

2 

0.110 

0.5  02 

0.467 

-0.035 

3 

0.150 

0.406 

0.452 

0.046 

4 

C.2C3 

0.4C2 

0.411 

C.009 

5 

C.250 

0.411 

0.391 

-0.020 

6 

0.275 

0.403 

0.390 

-0.013 

7 

0-.298 

0.398 

0.396 

-C.002 

8 

0-322 

0.411 

0.407 

-C.004 

9 

0,351 

0.421 

0.424 

0.003 

10 

0.372 

0.434 

0.439 

C.005 

11 

0.398 

0.443 

0.457 

0.014 

12 

0.427 

0.472 

G.477 

0.005 

13 

0.458 

0.501 

0.499 

-0.002 

14 

0.502 

0.554 

0.539 

-0.015 

15 

C.55  0 

0.622 

0.628 

C.  006 

Figure  16.   Velocity  vs  Thrust  Output 
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VELOCITY    VS    THRUST    6TH   ORDER    FIT 

COEFFICIENTS    OF    THE    POWER    SERIES    EXPANSION 

Ym  =  8(l)*e(2)*X  +  B(3)*X**2  +  fi(4)*X4*3+... 

B(  1)=  6.15428D-04  B(  2)=  1.7CC49D  CI  B(  3)=-2.03641D  02 
B(  4)=  1.C923CD  C3  B(  5-1  =-2.963230  03  B(  6)=  3.98231C  03 
B(    7)=-2. 099990    C3      B( 

ESTIMATES   OF    ERROR    FOR    THE    COEFFICIENTS 

ERRB(  1)=5.957D-C2  ERRB(  2)=7.585D  00  EFRB(  3J=1.348C  02 
ERRE(  4I=£.979D  02  ERR6(  5)=2.820D  03  ERRB(  6)=4.214D  03 
EPRB(    7)=2.415D    03       ERRB< 

SUN    SC,CEV    =    1..262D-04  F-RATIO   =    7.56CC-C1 

CHISQ=      4.978D-C3  DEG    CF    FREEDOM    =      8. 


DATA 

FIT 

CIFFERENCE 

I 

Mil 

Y(I) 

YY(  I) 

YY(II-Y(IJ 

1 

O.C 

0.0 

0.001 

0.001 

2 

0.110 

0.502 

0.486 

-0.014 

3 

0^150 

0.406 

0.434 

0.028 

4 

0.203 

C.4C2 

0.392 

-0.010 

5 

0.250 

0.411 

0.392 

-0.019 

6 

C..275 

C.4C3 

0.401 

^0.00  2 

7 

C.298 

0.398 

0.4C9 

C.011 

8 

0.322 

0.411 

0.418 

0.007 

9 

0.351 

0.421 

0.427 

0.006 

1C 

C.372 

0.434 

0.433 

-C.001 

11 

0.398 

0.443 

0.444 

0.001 

12 

0-.42  7 

0.472 

0.463 

-0.009 

13 

0.458 

0.501 

0.496 

-0.005 

14 

0.502 

0.554 

0.564 

0.010 

15 

G.550 

0.622 

0.619 

-0.003 

Figure    16o     (Continued) 
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VELOCITY  VS  THRUST  7TH  ORDER  FIT 

COEFFICIENTS  OF  THE  POWER  SERIES  EXPANSION 

Y(X)  =  B(l)J-e<2)*X«-fi(3)*X**2*fl(4)*X**3+... 

B(  1)=  1.138500-04  B(  2)  =  2.746530  CI  B<  2)=-4.453<7D  02 
B(  4)=  3.255C0D  03  Bt  5)=-1.26547D  04  B(  6)  =  2.71383C  04 
B(    7)=-3.02829D    C4       B<     8)=    1.373320    C4       B( 

ESTIMATES   OF    ERROR    FCR    THE    CCEFFICIENTS 

ERRB(  1)=3.674D-C2  ERRB(  21=1.1780  01  ERR8(  3)=2.634C  02 
ERRB(  4)=2.3040  C3  ERRB(  5)=1.017D  04  ERRB1  6J=2.408D  04 
ERRB(     7)=2.918D    04       ERRBt    8J=1.420D    04       ERRBC 

SUN    SQ    CEV    =    4-200D-05  F-RATIO    =    9.356D-C1 

CHISQ=       1.610C-C3  CE.G    CF    FREEDOM    =      7. 


-: DATA 

FIT 

CIFFERENCE 

I 

X(I) 

Yd) 

YYd) 

YYUJ-Ym 

1 

0.0 

0.0 

0.000 

0.000 

2 

0.11C 

0.502 

0.498 

-0.004 

3 

0.150 

0.406 

0.417 

0.011 

4 

0.203 

0.402 

0.392 

-0.010 

5 

C.250 

0.411 

0.4C4 

-0.007 

6 

0.275 

C.403 

0.408 

0.005 

7 

0-.2  98 

0.393 

0.410 

0.012 

8 

0-322 

0.411 

0.412 

0.001 

9 

Oi.351 

0.421 

0.417 

-0.004 

10 

0.272 

0.434 

0.425 

-0.009 

11 

0.398 

0.443 

0.444 

0.001 

12 

0.427 

0.472 

0.473 

0.001 

13 

0.458 

G.5C1 

0.508 

0.007 

14 

0..502 

0.554 

0.550 

-0.004 

15 

C55C 

0.622 

0.623 

0.001 

Figure  16.  (Continued) 
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VELOCITY    VS    THRUST    8TH    CRDER    FIT 

COEFFICIENTS    OF   THE    POWER    SERIES    EXPANSION 

Y(X)=B( l)*B(2)*X+B(3)*X*+2+B(4)*X**3+... 

B(  1)=  8.75C39D-06  B(  2)  =  4.58C88D  01  B(  3)=-9.48749D  02 
B(  4i=  8.817460  03  B(  5J=-4.496C7D  04  B(  6i=  1,342410  05 
B(    7)=-2.34717D    05       8(    8)=    2.222230    05       E(    9)=-8. 809140    04 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 

EPRB(  1)=1. 9730-02  ERRB(  2)=1.778D  01  ERRB(  3)=4.773D  02 
ERRBi  4)=5.188D  03  ERRB(  5)=2.976D  04  ERRB(  6)=9.795D  04 
EBRBl    7)=1.858D    C5      ERRB(     8)=1.890D    05      ERRB(    9)=7.S79D    04 

SUM    SQ    DEV    =    1.C38D-05  F-RATIO    =    1.219C    00 

CHISC=      3.735D-04  OEG    OF    FREEDOM    =       6. 


OAT  A 

FIT 

DIFFERENCE 

I 

•X(I) 

Yd  J 

YYU) 

YYUJ-YCI J 

1 

0.0 

0.0 

0.000 

0.000 

2 

0-11C 

0.502 

0.502 

-0.000 

3 

0.150 

0.406 

0.4C7 

0.001 

4 

0.203 

0.402 

0.401 

-0.001 

5 

C.250 

0.411 

0.4C8 

-0.003 

6 

0..275 

C.403 

0.4C6 

0.003 

7 

0.298 

0.393 

0.404 

0.006 

8 

0,322 

0.411 

0.407 

-0.004 

9 

0.351 

0.421 

0.418 

-0.003 

1C 

0-372 

0.434 

0.431 

-0.003 

11 

C..398 

0.443 

0.45G 

0.007 

12 

0.427 

0.472 

0.472 

0.000 

13 

0.458 

0.501 

0.498 

-0.003 

14 

0.502 

0.554 

0.555 

0.001 

15 

C.550 

0.622 

0.622 

-0.000 

Figure  16a  (Continued) 
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01 

e< 

3)=-l. 527640 

03 

04 

B< 

6)=  3.6^1510 

05 

06 

6( 

9)=-9.6G214D 

05 

VELCCITY  VS  THRUST  9TH  CRDER  FIT 

COEFFICIENTS  OF  THE  POWER  SERIES  EXPANSION 

Y(X)=B(1)*B(2)*X+B(3)*X**2+B(4)*X**3+... 

B(     1)=    9. 615780-08  B(  2)=    6.406660 

B(    4)=     1.64407C    C4  B(  5J=-9. 964990 

B(    7)=-8.55135C    05  B(  8)=    1.211070 

B(10J=    2.26633D    C5  B( 

ESTIMATES    CF    ERROR    FOR    THE    COEFFICIENTS 

EPRB(    l)=1.807D-02  ERRB(  2)=5.087D  01  ERRB(  3J-1.5S0C  03 

ERRB(     4)=2.C68D    04  ERRBt  5)=1.469D  C5  ERPEt  6)=6.263C  05 

ERRB(    7)=1.647D    06  ERRBt  8)=2. 6160  06  ERRB(  9J=2.3C3D  06 

ERRB(1C)=8.628D    05  ERRBt 

SUP    SC    CEV    =    7-2550-06  F-RATIO    =    1.425C-C1 

CHISQ=      2.610D-C4  DEG    CF    FREEDOM    =      5. 


DATA 

FIT 

CIFFERENCE 

I 

.xm 

Y(IJ 

YY(I) 

YYd)-Y(I) 

1 

o.c 

O.C 

0.000 

0.000 

2 

0-110 

0.502 

0.502 

0.000 

3 

0.150 

0.406 

0.405 

-0.001 

4 

C.203 

0.4C2 

0.404 

0.002 

5 

C.250 

0.411 

0.406 

-0.005 

6 

0.275 

G.403 

0.404 

0.001 

7 

0.293 

0.3S8 

0.404 

0.006 

8 

0.322 

0.411 

0.406 

-0.003 

9 

0.351 

0.421 

0.420 

-0.001 

10 

0.372 

0.434 

0.431 

^0.003 

11 

0.395 

0.443 

0.446 

0.005 

12 

0.427 

0.  472 

0.47C 

-0.002 

13 

0.458 

0.501 

0.5C1 

-C.000 

14 

0..  502 

0.554 

0.554 

0.000 

15 

C-550 

0.622 

0.622 

-0.000 

Figure    16c     (Continued) 
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00 

E( 

3)=  4.69988C 

02 

05 

B( 

6)=-9. 642410 

05 

C6 

B{ 

9i=  1.22215D 

07 

06 

B( 

VELCCITY  VS  THRUST  10TH  CRDEP  FIT 
COEFFICIENTS  OF  THE  POWER  SERIES  EXPANSION 
Y(X)=B(1)+E(2)*X+B(3)*X**2+B(4)*X<*3+... 

fi(  i)=-L 455160-07  B(  2)=  7.792420 

B(  4)=-l. 371370  04  B(  5) =  1.552C60 

B<  7)=  3,637360  C6  8(  8)  =  -8. 564500 

B<10)=-9. 921360  06  3(11)=  3.428700 

ESTIMATES  OF  ERROR  FCR  THE  COEFFICIENTS 

ERRB(  l)=1.777D-02  ERRB(  2)=2.077D 

ERRB(  4)=1.C99D  C5  ERRBC  5)=9.241D 

ERRE(  7)=i.617D  C7  ERRBl  8)=3.5110 

ERR8(1C)=3.671D  07  ERRB (  1 1  )  =  1  .2280 

SUP  SC  CEV  =  5-6150-06  F-RATIO  =  7.7840-02 

CHISQ=   1.96*0-04  CEG  CF  FREEDOM  =   4. 


02 

ERRB( 

3)=7.222D 

03 

05 

EPRB( 

6)=4.015C 

06 

C7 

ERRB( 

9) =  4. 762D 

07 

07 

ERRB( 

DATA 

FIT 

CIFFERENCE 

I 

xm 

Y(I) 

YY(I) 

YY(IJ-YU) 

1 

0.0 

0.0 

-0.000 

-0.000 

2 

0.11C 

G.5C2 

0.502 

C.  u 0 0 

3 

0.150 

0.406 

0.406 

-0.000 

4 

0.203 

0.402 

0.4G3 

0.001 

5 

C.250 

0.411 

0.40fc 

-0.003 

6 

0.275 

0.403 

0.404 

0.001 

7 

C.293 

0.398 

0.4C2 

0.004 

8 

0.322 

0.411 

0.407 

-0.004 

9 

0.351 

0.421 

0.421 

-0.000 

1C 

0*272 

C.434 

0.433 

-0.001 

11 

0.398 

0.443 

0.448 

0.005 

12 

0.427 

0.472 

0.46o 

-0.004 

13 

C.458 

C.5C1 

0.503 

0.002 

14 

0.502 

0.554 

0.554 

-0.000 

15 

C.550 

0.622 

0.622 

0.000 

Figure  16„  (Continued) 
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Figure    17.       Velocity    vs    Thrust    Curves. 
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Velocity   vs   Thrust   -   6th  Order  Fit 


Figure    17 


(Continued) 
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Velocity  vs   Thrust   -   7th  Order   Fit 


Figure    17.        (Continued) 
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Velocity  vs   Thrust   -   8th  Order   Fit 


Figure    17.        (Continued) 
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Velocity  vs   Thrust   -   9th  Order   Fit 


Figure    17 


(Continued) 
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Figure    17.        (Continued) 
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SYM 


SYM   0 


SYM   1 


SYM   2 


APPENDIX    C 
FLOW    CHART    FOR    DATA    REDUCTION    PROGRAM 


1.  depress  JUMP,    0,    0 

2.  load   card    "DATA  REDUCTION   ROUTINE" 

3.  place   sense   switch  up  for  point  plot 

4.  depress  PEN  UP   switch  on  plotter 


z 


1.  depress  JUMP,  1,  3 

2.  place  RUN/STEP/LOAD  switch  to  LOAD 

3.  enter  code  for  y  axis  channel* 

4.  depress  LIST  PROG  4  times 

5.  enter  code  for  x  axis  channel* 

6.  put  RUN/STEP/LOAD  switch  to  RUN 


CHANNEL   CODE 

1 340 

2 140 

3 240 

4 040 

5 300 

6.,... 100 

7 200 

8 000 


1.  depress  JUMP,  SYM,  0,  RESUME 

2.  set  0.00  (button  left  corner) on  grid  with  ZERO  controls 

3.  depress  RESUME 

4.  set  10.00  (top  right  corner) on  grid  with  SCALE  controls 

5.  enter  X  grid  coordinate  of  origin,  depress  RESUME 

6.  enter  Y  grid  coordinate  of  origin,  depress  RESUME 

7.  pen  moves  to  origin,  depress  REMOTE  switch  on  plotter 


NO 


depress   JUMP, 
SYM,    5,    RESUME 


YES 


© 


1. 

2. 
3. 

4. 

enter  X  max,  RESUME 
enter  X  min,  RESUME 
enter  Y  max,  RESUME 
enter  Y  min,  RESUME 

' 

' 

1. 
2. 
3. 

enter  X  scale  (units/0.1  in),  RESUME 
enter  Y  scale  (units/0.1  in),  RESUME 
read  length  of  x  and  y  axes 

--f^      YES   v 

depress  JUMP 

CVM     9     DPCTTMC 

NO 
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SYM  3 


1. 

enter  interval  between  graduation  of  X,  RESUME 

2. 

enter  interval  between  graduation  of  Y,  RESUME 

axes  are  plotted 

E 


depress  JUMP,  SYM,  5,  RESUME 


SYM  5 


1.  enter  x  scale  (in/volt)  for  plot,  RESUME 

2.  enter  y  scale  (in/volt)  for  plot,  RESUME 

3.  depress  FLAG  if  printout  is  desired 

4.  enter  X  scale  (units/5  volts)  for  printout,  RESUME 

5.  enter  Y  scale  (units/5  volts)  fcr  printout,  RESUME 


next  point  is  plotted 
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APPENDIX    D 
DATA    REDUCTION    PROGRAM    LISTING 


Statement  Number 


Octal 

Decimal 

Code 

1 

0 

0000 

266 

1 

1 

1  60 

2 

2 

220 

-3 

3 

0  1  6 

4 

4 

266 

5 

5 

1  6  1 

6 

6 

220 

7 

7 

0  1  6 

10 

8 

266 

1 

9 

1  62 

2 

00  1  0 

220 

3 

1 

0  1  6 

4 

2 

266 

5 

3 

1  63 

6 

4 

1  60 

7 

5 

3  1  4 

20 

6 

260 

1 

7 

025 

2 

8 

3  1  4 

3 

9 

220 

4 

0020 

034 

5 

1 

302 

6 

2 

247 

7 

3 

025 

30 

4 

16  1 

1 

5 

3  1  4 

2 

6 

24  1 

3 

7 

05  1 

4 

8 

360 

5 

9 

326 

6 

0030 

3  1  5 

7 

1 

273 

Symbol 


a* 


n 


<£ 


H 


Explanation  (if  needed) 


load  IX  register 
with  160 
jump  to  octal 
address  16 

load  IX  with 

161 

jump  to  octal 

address  16 

load  IX  with 

162 

jump  to  octal 

address  16 

load  IX  with 

163  J 

store  IX  into  IRO  register 

OPCB  (output  control  byte) 

if  IL    of  IX  ?   0  jump 

to  octal  address  25 

OPCB 

jump  to  octal 

address  34 

RIIJS  (read  input  status) 

if  H8  of  IX  =  1,  jump  to 

octal  address  25 

recall  IRO  register  to  IX 

OPCB 

if  L2  of  IX  =  1,  jump  to 

octal  address  51 

clear  IX 

transfer  digit  of  E  to  IL 

O^DE    (sign  transfer) 

load  DH  with  111  111 


to  transfer 

data 

to  plotter 

to  transfer 
Y  data 


pen  up 


pen  down 
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40 

0032 

277 

load  DL  with  Oil  111  111 

1 

3 

322 

-  load  IX  with  (DC)* 

2 

A 

3  1  5 

OPDB  (transfer  exponents) 

3 

5 

336 

decrement  DC  by  1  step 

4 

6 

322 

load  IX  with  (DC)* 

5 

7 

3  1  5 

OPDB  (transfer  first  2  digits) 

6 

8 

336 

decrement  DC  by  1  step 

7 

9 

322 

load  IX  with  (DC) * 

50 

0040 

3  1  5 

OPDB  (transfer  next  2  digits) 

1 

1 

057 

RESUME  (end  plot  subroutine) 

2 

2 

066 

"^  symbolic  52:  subroutine  to 

3 

3 

052 

F 

J  connect  binary  to  decimal 

4 

4 

1  1  4 

* 

clear  data  register 

5 

5 

000 

0 

0 

6 

6 

063 

a 

clear  entry  register 

7 

7 

377 

no  operation 

60 

8 

240 

jump  if  Ll  =  1  to 

1 

9 

1  20 

* 

octal  addres  120 

2 

0050 

24  1 

jump  if  L2  =  1  to 

3 

1 

1  27 

Br 

octal  address  127 

4 

2 

242 

jump  if  L.4  =  1  to 

5 

3 

1  35 

octal  address  13  5 

6 

4 

243 

jump  if  LS  =  1  to 

7 

5 

1  43 

octal  address  143 ' 

70 

6 

244 

jump  if  Hi  =  1  to 

1 

7 

1  5  1 

octal  address  151 

2 

8 

245 

jump  if  H2  =  1  to 

3 

9 

1  56 

octal  address  155 

4 

0060 

246 

jump  if  H4  =  1  to 

5 

1 

1  63 

octal  address  163 

6 

2 

24  7 

jump  if  HS  =  1  to 

7 

3 

1  70 

octal  address  170 

100 

4 

1  1  1 

t 

1 

5 

000 

0 

^ 

2 

6 

024 

+ 

to  reconvert 

3 

7 

002 

2 

input  data 

4 

8 

005 

5 

E  register  now 

5 

9 

006 

6 

r  contains  input 

6 

0070 

020 

s 

J   data  valine 

7 

1 

057 

1  RESUME 

82 


110 

007  2 

1 

3 

2 

4 

3 

5 

4 

6 

5 

7 

6 

8 

7 

9 

120 

0080 

1 

1 

2 

2 

-3 

3 

4 

4 

5 

5 

6 

6 

7 

7 

130 

8 

1 

9 

2 

0090 

3 

1 

4 

2 

5 

3 

6 

4 

7 

5 

140 

6 

1 

7 

2 

8 

3 

9 

4 

0  100 

5 

1 

6 

2 

7 

3 

150 

4 

1 

5 

2 

6 

3 

7 

4 

8 

5 

9 

6 

0  110 

7 

1 

377 
377 
377 
377 
377 
377 
377 
377 
00  1 
002 

0  1  0 

1  1  3 
000 
220 
062 
006 
004 
1  1  3 
000 
220 
064 
003 

00  2 

1  1  3 
000 
220 
066 

00  1 
006 

1  1  3 
000 
220 
070 

0  1  0 

1  1  3 
000 
220 
07  2 
004 
1  1  3 


J 
2 


e 


R° 


-\ 


no  operations 


"^add   128   to  main  data   register    (MDR) 
0   if  most   significant 
but    (Ll)    of   IX  register 
is   on    (Ll   =   1) 

jump  to   octal 

address  62 
'  add  64  to  MDR  0 
if  L2  =  1 


jump  to  octal 
address  64 


> 


add  32  to  MDR  0 
if  L4  =  1 


jump  to  octal 

address  66 
~> 

,add  16  to  MDR  0 
if  L8  =  1 


jump  to  octal 

address  70 
~")add  8  to  IIDR  0 
Jif  III  =  1 

jump  to  octal 

address  72 


3 


add  4  to  MDR  0 
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160 

0  112 

000 

1 

3 

220 

2 

4 

074 

3 

5 

002 

4 

6 

1  1  3 

5 

7 

000 

6 

8 

220 

7 

9 

076 

170 

0  1  20 

00  1 

1 

1 

1  1  3 

2 

2 

000 

3 

3 

220 

4 

4 

1  00 

5 

5 

377 

6 

6 

377 

7 

7 

377 

200 

8 

377 

1 

9 

266 

2 

0  1  30 

000 

3 

1 

3  1  4 

4 

2 

305 

5 

3 

163 

6 

4 

266 

7 

5 

000 

210 

6 

3  1  4 

1 

7 

305 

2 

8 

057 

3 

9 

066 

4 

0  1  40 

075 

5 

1 

200 

6 

2 

0  1  4 

7 

3 

200 

220 

4 

0  1  0 

1 

5 

056 

2 

6 

057 

3 

7 

066 

4 

8 

070 

5 

9 

023 

6 

0  150 

1  1  1 

7 

1 

002 

u* 


if  H2  =  1 
jump  to  octal 

address  74 

add  2  to  MDR  0 
'  if  H4  =  1 

jump  to  octal 
address  76 

add  1  to  MDR  0 

if  H8  =  1 

jump  to  octal 
address  100 

•no  operations 
next  is  channel  select  and 

data  transfer  code  for  A/D 

converter 
load  IX  with 

(code  for  Y  data  channel) 

OFCB)get  y  value  in  IX 

IPDB) 

store.  IX  in  IR1  register 

load  IX  with 

(code  for  x  data  channel) 

0PCB)get  x  value  in  IX 

IPDB) 

RESUME 

}  symbolic  75  point 
plot  subroutine 

branch  to  octal 
address  14 
branch  to  octal 

address  10 
halt 

RESUME 

}  symbolic  sin  x  data  scaling 
and  transfer  subroutine 

•EX  MDR  2  for  x  scale 
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230 

0  152 

02  1 

♦ 

x 

1 

3 

1  1  1 

t 

>  add  x  origin 

2 

4 

0  1  0 

•  9 

3 

5 

020 

s 

J 

4 

6 

220 

"\  jump  to  x  data  entry  of 

5 

7 

000 

0 

fplot  routine 

6 

e 

066 

-\  symbolic  arc  sine  Y  data 
_J  and  transfer  subroutine 

scaling 

7 

9 

07  1 

5 

240 

0  160 

02  3 

X 

^ 

1 

1 

1  1  1 

t 

SEX  MDR  3  for  Y  scale 

2 

2 

003 

J 

J 

3 

3 

02  1 

♦ 

4 

4 

111 

t 

y.   add  Y  oriain 

5 

5 

0  1  1 

9 

\ 

6 

6 

020 

s 

J 

7 

7 

220 

~)    jump  to  Y  data  entry 

250 

8 

004 

4 

J   plot  routine 

1 

9 

066 

~N  symbolic  5  scale  factors 

2 

0  170 

00  5 

S 

[  storing  subroutine 

3 

1 

065 

advance  paper 

4 

2 

1  66 

clear  flag 

5 

3 

1  1  7 

set  decimal  point  to 

6 

4 

003 

J 

3  places 

7 

5 

200 

n  jump  to  pen  up  entry 

' 

260 

6 

0  1  0 

# 

1   of  n)lot  routine 

1 

7 

056 

halt 

. 

2 

8 

1  10 

♦ 

"] 

3 

9 

002 

2 

>  enter  and  print  x  scale 

4 

0  180 

060 

J   for  plot 

5 

1 

056 

halt 

6 

2 

1  1  0 

* 

enter  and  print  x  scale 

7 

3 

003 

J 

V  for  printout 

270 

4 

060 

J 

1 

5 

056 

halt 

2 

6 

1  1  0 

1 

I  enter  and  print  Y  scale 

3 

7 

004 

4 

>  for  plot 

4 

8 

060 

J 

5 

9 

056 

halt 

6 

0  190 

1  1  0 

* 

V  enter  and  print  Y  scale 

7 

1 

005 

S 

J 

85 


300 

0  192 

1 

3 

2 

4 

3 

5 

4 

6 

5 

7 

6 

8 

7 

s9 

310 

0200 

1 

1 

2 

2 

3 

3 

4 

4 

5 

5 

6 

6 

7 

7 

320 

a 

1 

9 

2 

02  10 

3 

1 

4 

2 

5 

3 

6 

4 

7 

5 

330 

6 

1 

7 

2 

8 

3 

9 

4 

02  20 

5 

1 

6 

2 

7 

3 

340 

4 

1 

5 

2 

.6 

3 

7 

4 

8 

5 

9 

6 

0230 

7 

1 

060 
06  5 
066 
006 
200 
20  1 
200 
054 
1  2  7 
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TABLE  VII 
DATA  REDUCTION  PROGRAM  SCRATCH  PAD  REGISTER  USAGE 
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USAGE 
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AXES  DRAWING  SECTION 

PLOTTING  SECTION 

X  maximum 

decimal  conversion 
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X  minimum 

temporary  storage 
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X  scale  (units/ 0.1  in) 

X  scale  (V/in) 

3 

Y  scale  (units/0.1  in) 

Y  scale  (V/in) 
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Y  maximum 

X  scale  (units/V) 
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Y  minimum 

Y  scale  (units/V) 
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X  interval 

not  used 
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Y  interval 

not  used 
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grid  origin  X 

grid  origin  X 
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grid  origin  Y 

grid  origin  Y 
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